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ABSTRACT 

In  this  paper  we  develop  time-domain  state  space 
models  for  lossless  layered  media  which  are  des- 
cribed by  the  wave  equation  and  boundary  conditions. 
Our  models  are  for  non-equal  one-way  travel  times; 
hence, they  are  more  general  than  existing  models  of 
layered  madia  which  are  usually  for  layers  of  equal 
one-way  travel  times.  Fuil  state  models, which  In- 
volve 2K  states  for  a K-layar  media  system, as  well 
as  half-stata  models, which  involve  only  K states 
are  developed  and  related.  Certain  transfer  func- 
tions,which  appear  In  the  geophysics  literature  in 
connection  with  models  of  layered  media  with  equal 
travel  times, are  general  ized  to  the  situation  of  non- 
equal travel  times. Our  state  space  models  represents 
new  class  of  squat  Ions,  causa  I functional  aquations , 
some  of  whose  properties  and  approaches  to  simulation 
are  discussed. 

I.  IMTROOUCTIOW 

We  are  interested  in  lossiess  layered  madia  which 
are  described  by  the  wave  equation  and  boundary  con- 
ditions. Specific  applications  of  such  media  are: 

(1)  horizontaily  stratified  nonabsorptiva  earthwith 
vertically  traveling  plane  compress i ona I waves, and 

(2)  interconnection  of  lossless, not  necessarily 
matched, transmission  lines.  Other  applications  can 
be  found  in  acoustics, optical  thin  coatings, EM  prob- 
lems,etc.  Because  our  interest  is  presently  in  the 
seismic  area, our  discussions  will  be  in  the  context 
of  such  models. 

A system  of  K layered  media  is  depicted  in  Fig. 1.  We 
adopt  the  convention  of  tailing  the  layer  below  lay- 
er K the  basement.  The  basement  is  assumed  to  act 
like  an  energy  sink;  l.e.,no  energy  is  returned  from 
the  bascaent  into  the  K layers.  Each  layer  is  char- 
acterized by  its  one  way  travel  time.T;,  velocity ,Vt, 
and  normal  Incidenca  reflection  coefficient  ri  (iw 
1,2,...,K).  Additionally  Interfaca-O  denotes  the 
surface  and  is  characterized  by  reflection  coeffi- 
cient rg.  In  Flg.t,ffl(t)  and  y(t)  denota  the  Input 
(e.g. ,sal smic  source  signature  from  dynamlte,ai rgun, 
etc.)  to  the  layered  media  system  which  is  applied 
at  intarfaca-0,  end  the  output  (i.e., ideal  seismo- 
gram) of  the  system  which  is  observed  at  the  surface, 
respectively.* 

We  shall  prasent  state  space  models  for  the  Flg.1 
system.  These  models, as  we  shall  describe  more  fully 

*ln  a marine  environment, layer  I can  be  taken  to  be 
water;  but, In  that  case  m(t)  is  applied  and  y(t)  Is 
observed  In  the  water  layer.  It  is  relatively 
straightforward  to  extend  the  results  of  this  paper 
to  that  ease. 


below, are  quite  different  from  those  which  have 
appeared  in  the  Geophysics  iiteraturq  (Refs. I -4, 
for  example).  One  big  difference  Is  that  ourmodels 
are  for  non-equal  one-way  travel  times. 

An  Important  use  of  a modal  of  a K-layer  madia  sys- 
tem Is  to  generate  synthetic  seismograms;  i.e., to 
generate  y(t)  for  a given  m(t).  This  synthetic 
data  can  then  be  used  for  preliminary  testing  and 
evaluation  of  signal  processing  techniques  (e.g., 
deconvolution).  These  models  may  also  be  useful 
for  identifying  important  parameters, such  as  re- 
flection coefficients  and  one-way  travel  times.** 

As  in  Refs. 2 and  3, we  shall  find  it  convenient  to 
draw  ray  diagrams  with  time  dlsplacaewnt  along  the 
horizontal  axis,  so  that  the  rays  appear  to  be  at 
non-normal  incidence  and  so  do  not  overlap  one  an- 
other. Figure  2 depicts  primary  and  multiple  re- 
flections for  a 2-layer  madia  system, and  illus- 
trates the  very  complicated  internal  behavior  of 
even  a 2-layer  system.  Our  state  space  models  can 
not  only  be  used  to  compute  y(t) .but  can  also  be 
used  to  compute  the  internal  behavior  of  a layered 
madia  system.  They  are  based  on  ray  theory, which 
gives  exact  results  for  loss  less .horizintal ly  stra- 
ti f i ed  med i a . 

The  starting  point  for  our  developments  is  Che  Fig. 
3 ray  diagram.  Symbols  u^  and  d^  denote  the  upgo- 
Ing  and  downgoing  waves  in  Che  kch  layer, respec- 
tively; and.we  adopt  the  convention  chat  waves  at 
the  top  of  a layer  occur  at  present  time.t.  That 
each  layer  is  characterized  by  two  signals  travel- 
ing In  opposite  directions  follows  directly  from 
Che  solution  of  a lossless  wave  equation.  Geophys- 
icists will  recognize  chat  Fig. 3 Is  also  the  start- 
ing point  for  the  models  which  appear  in  Refs. I 
through  4, for  example.  We  shall  return  to  this 
point  shortly. 

As  stated  by  Robinson  (Ref. 3),  "the  solution  of  the 
wave  equation  at  each  Interface  leads  to  the  defi- 
nition of  a raflectlon  coefficient  rj  associated 
with  that  interface.  ...the  reflection  coefficient 
rj.which  must  satisfy  |rj|<l,has  these  properties. 
A downqolnq  wave  of  amplitude  A In  layer  J.upon 
striking  interface  j.  Is  both  reflected  and  trans- 
mitted. The  reflected  portion  Is  an  upgoing  waveof 
amplitude  rjA  In  layer  j,so  rj  represents  the  ra- 
flectlon coefficient.  The  transmitted  portion  i s a 


**Tha  reader  interested  in  elements  of  the  seismic 
prospecting  method  and  the  seismic  reflection 
technique  should  see  Ref.S.Chs.l  and  3, and  Refs. 
6,  7,  and  3. 
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downgoing  w«v«  of  amplitud*  (1*r,)A  in  layar  j+I,»o 
)+Tj  reprasents  eha  eranamission'eoaff Iciant.  An 
upqoinq  wava  of  aaiplitudaB  in  layer  J'^1  is  both  re- 
rlacted  and  cransmicted  when  ic  strikes  interface  j. 
The  reflected  portion  is  a downgoing  wave  of  ampli- 
tude Tj S in  layar  J+l,  and  the  transmitted  portion 
Is  an  upgoing  wave  of  amplitude  (l•rJ)B.  Hence -r^ 
and  (iTj)  represent,  respectively. the  reflection 
coefficient  and  transmission  coefficient  /or  the 
upgoing  wave.  These  properties  are  summarized  in 
Table  I (2)." 

Table  1.  Reflected  and  Transmitted  Portions 


Downgoing  wave 
A 

in  layer  J 

Upgoing  wave 
S 

in  layer  j+1 


Waveform  u^(t'«T^)  (Figure  3)  Is  made  up  of  two 
parts, namely  the  part  due  to  the  reflected  portion 
of  djU'tk)  and  the  part  dua  to  the  transmitted 
portion  of  u^^(t).  It  satisfies  the  equation 


Reflected 

Port  ion 

Transmi tted 
Port i on 

Upgoing  wava 
r^A 

in  layer  j 

Oowngoing  wava 
(Hr  )A 

In  layer  j-)-l 

Oowngoing  wave 
-r^B 

In  layer  j-*1 

Upgoing  wava 
(l-r  )B 

In  layer  j 

In  a similar  manner, waveform  d«.i (t)  satisfies  the 
equation 

■ (l-hr.  )d'(t-T,  ) - r.u._^,(t).  (2) 


We  refer  to  Eqs.  (1)  and  (2)  as  the  Interface  equa- 
tions. These  equations  are  the  starting  point  ^or 
transfer  function  models, which  are  very  popular  In 
the  Geophysics  literature  (Ref .2) .and, they  are  also 
our  starting  point  for  the  development  of  time- 
domain  state  space  models  (sea  Fig, A). 

In  the  Geophysics  literature,  the  assumption  of 
eoual  one-way  travel  times, T|, is  usually  made.  Lay- 
ers of  different  travel  times  are  built  up  by  In- 
serting layers  whose  reflection  coefficients  are 
zero.  Transfer  functions  tdiich  relate  the  upgoing 
and  downgoing  signals  from  one  layer  to  the  next 
have  been  obtained.  Because  of  the  equal  one-way 
travel  time  assumption, soma  vary  clever  recursive 
equations  have  been  developed  (Rafs.3,‘*,and  9,  for 
example)  which  simplify  the  calculations  of  the 
transfer  functions  relating  layers.  These  transfer 
functions  are  often  presented  in  terms  of  z-trans- 
forms, which  is  again  a consaquenca  of  the  equal 
travel  time  assumption.  From  the  transfer  functions 
reiating  layers, it  is  possible  to  obtain  Che  so- 
called  refleetJ[ontrBnsfer_funct]an,  Y(s)/M(s)  [or 
Y(z)/N(z} ] .which  is  used  to  generate  synthetic  seis- 
mogram data  (Rafs.3,^,and9,  for  example).  Addition- 
al ly, Robinson  (Ref. 10)  has  presented  a vary  simple 
technique  for  computing  Y(z)/H(z)  in  a layer-recur- 
sive manner.  Inverse  problems  of  reconstructing  , the 
reflection  coefficients  from  knowledge  of  y(t)  and 
m(t)  have  been  extensively  studied.  Their  solutions 
(Refs. 3, i*. and  II , for  example)  are  again  quite 
strongly  dependent  on  the  equal  travel  time  assump- 
tion,and  usa  z-transform  relationships. 

As  Is  well  known, there  is  a vast  literature  associ- 


ated with  systems  which  are  described  by  time- 
domain  state  space  models.  Host  recent  results  in 
estimation  and  identification  theory, for  example, 
require  a state  space  model.  In  fact, our  ultimate 
objective  is  to  apply  those  theories  to  the  layered 
madia  problem;  but, to  do  so, of  course  requires 
>T*te  space  models.  One  might  argue  that  it  should 
be  possible  to  go  directly  from  the  transfer  func- 
tions,already  developed, to  equivalent  state  space 
models.  In  most  cases, this  is  not  practical  since 
closed-form  expressions  for  the  reflection  trans- 
fer function  are  not  available;  that  transfer  func- 
tion must  be  computed  from  a set  of  equations  which 
are  solved  in  a recirsiva  manner.  Addi t ional ly , the 
transfer  function  results  which  appear  in  tha  Geo- 
physics literature  are  limited  by  the  equal  travel- 
time assumption. 

In  this  paper ,we  davelop  (sea  Fig. A)  a variety  of 
state  models  for  tha  FIg.1  system.  Full -state  mo- 
dels.which  are  of  dimension  2K  (i.a.,two  states 
per  layer) .are  described  In  Sections  11  and  V.A. 
Half-state  models, which  are  of  dimension  K (I.e., 
one  state  par  layer)  .are  described  In  Sections  IV 
and  V.B.  A number  of  useful  transfer  function  re- 
lationships are  described  In  Section  III.  Those  re- 
lationships are  obtained  directly  from  state  aqua- 
tions .and, not  only  serve  to  connect  our  results 
with  tha  transfer  functions  found  in  the  Geophys- 
ics I i terature.but  also  represent  generalizations 
of  those  results  to  the  non-equal  travel  time  case. 
Our  state  equations  are  continuous  time  equations 
with  multiple  t ima-delays .and.are  referred  to  at 
causal  functional  equations.  Some  preliminary  dis- 
cussions  on  simulation  o^  such  equations  Is  given 
In  Section  VI.  Some  new  results  which  have  been 
obtained  using  our  state  space  models,  as  well  as 
directions  for  future  research, are  given  in  Sec- 
tion VII. 

I I . A STATE  EQUATION  HODEL 


A state  equation  modai  for  our  K layer  media  sys- 
tem Is  obtained  directly  from  Eqs.(1)  and  (2), 
which  are  applicable  for  kw| ,2, . . . ,K-I , and  com- 
parable equations  at  the  surface  and  Kth  interface. 
At  tha  surface  (F!g.Sa),we  obtain 

y(t)  ■ rgm(t)  + (I-rg)uj(t)  (3) 

dj(t)  « (1+rg)  m(t)  - rgU^(t);  (A) 

and, at  the  Kth  interface,  we  assume*  that  u,^^  (t)w 
0,tP  obtain  (Fig. 5b) 


u^(t*T^) 


Signal  y(t)  In  Eq.(3)  is  tha  measurable  systemout- 
put.  Signal  d|^«i(t)  is  also  a system  output;  but, 
since  It  cannot  be  measured, we  shall  Ignore  It  in 
following  analyses. 

It  is  convenient  to  group  Eqs . ( 1) , (2) , (A) ,and  (5) 

In  a layer  ordering. as  follcxs: 


*This  assumption  Is  a boundary  condition  which  Is 
compatible  with  what  is  meant  by  the  "basement." 
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dj  (l)— r^Uj  (t)+(t+rg)m(t) 
u,(t*T,)-rjdJ(t-T,)*(t-r^)u2(?) 
dj(t)-(..rj.,)dj.,(t-Tj.,)-r..,Uj(t))  

d j;, ( t ) -( 1 I ) d 1^. , ( j I U|^(  t ) 

This  systsm  of  2K  equations  Is  not  In  a useful 
state  equation  format,  yet,  since  signals  In  Its 
left-hand  side  occur  at  t and  delayed  t lnMS,and s Ig- 
nals  on  the  right-hand  side  occur  at  t,  t-Tj.^and 
C-Tj . In  order  to  put  Eq.(7)  Into  a useful  state 
equation  forniat,  let 


dj(t)  S <lJ(t-Tj) 


for  all  j>l,2,...,K.  Observe, from  Fig. 3,  that  the 
dawngoing  states  d)(t)  occur  at  the  bottom  of  a 
layer.  Equation  (7)  becomes 

dj  (tyTj)— r^Uj  {t)+'{l+rg)m{t) 
u, (t*Tj)-r|dj (t)+(l-rj)uj(t) 
d.(t^.)-(i.r..,)d^.,(t)-r..,u^(t))  

Uj(t'*Tj.)-rjdj(t)f(l-r^.)uj^,  (t)  ) 

d^(t4^^)-(l^r^.,)dK-,(t)-r^-,u^(t) 

9y  means  of  transformation  (8)  each  pair  of  aqua- 
tions in  (7)  now  only  involves  two  time  points, 
t>T.  and  t.  Equations  (9)  and  (3)  together  repre- 
sent the  state  aquation  model  for  the  output  y(t). 
This  model  is  referred  to  as  the  layer-ordered 
(L-0)  full-state  model  In  the  sequel . 

Equations  (9)  and  (3)  can  be  expressed  in  more  com- 
pact notation  by  introducing  tha  following  2K x 2K 
matrix  operator*: 

^ - diag(z^,z,  ,22,i2,...,Z|^,*n),  (10) 

where  Zi  is  a scalar  operator  used  to  denote  a t, 
sec.  time  delay  (l.e.,Z|f(t)  • f(t-Tt)).  Let 

X(t)«col  (u^  (t)  ,dj  (t)  ,Uj(t)  ,dj(t) U|^(t)  ,d|^(t))  j 

(11) 

then,Eqs.(9)  and  (3)  can  bn  written, as 

y(t)  - c'i(t)  ♦ rgm(t)  (13) 

where  the  explicit  structures  of  A,  and  c,  can  ba 
deduced  directly  from  tha  former  aquatlons.~Because 
we  do  not  need  this  Information  at  this  point, and 
because  other  versions  of  Eqs.(12)  and  (13), which 
we  discuss  in  Section  V-A,have  mere  easily  remem- 

»tliis  idea  was  first  suggested  to  us  by  Hr.HIchael 
Steinberger,a  graduate  student  In  the  Electrical 
Engineering  Department , at  the  University  of 
Southern  California. 


bared  A,  and  £matricas,  we  do  not  give  explicit 
A,  £,  and  c structures  here  for  Eqs.(12)  and  (13). 
From  Eqs.(T2)  and  (13).  we  see  that 

i(t)-(?'’-A)*'^(t)  - (l-;A)'';bm(t)  (lit) 
and  ~ 

y(t)-[c‘(;''-A)‘Vrg]m(t)-tc'(l-;A)‘’pb+rg]m(t). 

(IS) 

Thesa  equations  provide  us  (conceptually,  at  least) 
with  the  solution  of  tha  state  equation  and  with 
the  output  as  a function  of  the  Input.  The  trans- 
fer function  of  the  K layer  media  system  is  ob- 
tained  directly  from  Eq.(l5).  as 

^ -£•(;■' -A) - c-(l-^)*’;Vrg  (16) 
where  p Is  obtained  from  p by  setting 


At  this  point, soma  comments  on  Che  nature  of  stats 
equation  (12)  are  in  order.  In  the  special  case 
when  i T,p  “ z I, where  z denotes  the 

T sac  time  delay  and  1 is  the  2Kx2K  identity  ma- 
trix; and,Eq.(l2)  can  be  written, as 

X(t'hT)  • A^(t)  * bm(t)  . (18) 

This  equation  can  be  reduced  to  a vector  finite- 
difference  equation  by  choosing  twkT, and, when  m(t) 
only  has  values  at  c^kT.  Then, all  of  the  usual 
techniques  associated  with  such  equations  can  be 
used  CO  analyze  our  K-layer  media  system.  We  donot 
choose  CO  follow  this  uniform  travel  time/sampled 
data  path, because  these  assumptions  seem  coo  re- 
strictive. 

State  Equation  (12)  is  a dynamical  equation  with 
multiple  time  delays.  It  is  not  a differential 
equation, nor  Is  it  a finite-difference  equation. 

We  shall  refer  to  it  as  a causal  functional  equa- 
tion. It  is  linear  and  t i me- 1 nvar i ant ,and , as  I s the 
case  with  delay-time  systems ,requi res  initial  val- 
ue Information  over  initial  intervals  of  time. 
Equation  (16)  suggests  a straightforward  way  to 
compute  y(t)  for  an  arbitrary  m(t).  First  compute 
the  system's  impulse  response, H(s) .where, obviously 

H(s)  • £'(l-^)*’5b+rg 

then, convolve  h(c)  with  m(c)  to  obtain  y(c).  It  is 
Interesting  to  note  that  h(t)  is  a sequence  of  im- 
pulse functions. since  t he  right-hand  side  of  Eq. 

( 19)  is  an  infinite  series  each  of  whose  terms 
looks  like  ae'*®,and„^"Kae**®}  • a6(t-B). 

Since  our  K layer  media  system  is  one  with  time 
delays, Its  state  space  is  Inf ini te-dimens ional ;but 
interestingly  enough, only  a finite  number  of  states 
(l.a.,2K)  are  needed  to  describe  the  trajectories 

**The  rlght-fiand  side  of  Eq.(t7)  is  the  Laplace 
transform  of  a delay  effect.  Associating  z,  with 
a delay  is  common  In  the  Seophysics'  and  time- 
series'  literatures.  In  the  control's  literature, 
on  the  other  hand.z,  is  usually  associated  with 

st; 

an  advance  in  which  case  z, ~e  . Clearly, our  z. 

Is  simply  the  Inverse  of  z^  in  the  control's  lit- 
erature. 
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In  chat  scat#  space. 

Tha  important  systam  thaoratic  eoncapcs  of  obsarva* 
bi 1 1 ty .control labi I i ty .and  idant i f iabi 1 1 Cy  can  ba 
dafinad  in  a variaty  of  ways. as  for  di ffarant lal * 
dalay  aquations, and  xill  ba  discussad  aisawhara. 

Tha  L*0  full-scata  model  was  first  prasancad  by 
Nahi  and  Hendal  In  Raf.I2. 


111.  SOHE  TRANSFER  FUNCTIONS 

Whila  Eg. (16)  Is  tha  transfer  function  for  our  K* 
layer  madia  systam. it  is  not  at  all  a useful  form 
for  computing  that  quantity.  In  this  saction  wa  pre- 
sent two  alcarnativas  to  Eq.(l6),both  of  which  are 
recursive  in  nature  and  are  of  interest  in  thair 
Okm  right. 

A.  Laver  Transfer  Functions 

Lac  £,(s)  denote  a 2xi  Laplace  transformed  vector, 
dafinad  for  tha  kth  layer,  as 

X^^(s)  » col  [U|^(s)  .0|^(s)]  . , (20) 

In  Appendix  A,  wa  show  that 

X.(s)  - a [W  '(s)G,(s)]M(s)  (21) 

i-k  ' 
and 

Jlfy-  (l-rg)(1.0)w'’(s)6,(s)v.rg  (22) 


where  Uj(s)  is  a 2X2  matrix  which  is  solved  in  a 
backwards  racursiva  manner  from  tha  following  algor- 
i thm: 


^^{s)  - ^^(s) 

Wj(s)  - F|(s)-H,(s)w‘],C,^,(s) 
I - K-1,IC-2 1 

where 


(23) 


and 


(24) 

(25) 

(26) 

(27) 


In  practice. It  Is  not  necessary  to  compute  the  four 
alamants  of  W*'^(s),duo  to  tha  sparse  nature  of 
H,  (s)  and  G„i(s).  Only  tha  1-2  element  of 
is  non-zero,  and  that  alemant  depends  only  on  tha 
1-2  element  of  W^iiCs). 

During  tha  development  of  Eq.(21)we  obtain  l^j(s)  as 
a function  of  l<4-i(s): 


X.(s)  - w‘'(s)Gj(s)Xj.,(s)  , j-2,3 K 

-1 

X^ (s)  • (s)Gj (s)H(s) 

Matrix  uT‘(s)G.(s)  is  tha  layer  transfer  function 
matrix  which  relates  layers  J-l  and  j.  Equation 
(28)  Ts  similar  to  tha  racursiva  algorithms  which 
appear  in  tha  Geophysics  literature  (Refs. 3, 4 and 
9, for  example) .except  that  it  is  a generalization 
of  those  algorithms  to  non-equal  travel  times. 


We  conclude  this  paragraph  with  an  example  which 
illustrates  tha  calculation  of  Y(s)/M(s)  for  Kw2. 
In  that  case. 


(29b) 


W,-fi-HiW-'G2 


-(r£22^ 

’**•1^2 


(29c) 


W,  (1,2) 


and 


’"•'•i'’2V''0^1*^'-0"2^M 


mTs 


rg*(l-r^)z,w;'(l.2) 


(29d) 


2 2 2 2 
Vo"!  Wl*l*'’2*l*2 
1 r~^  1 TT 

'*'’lWo'’l*l*V2*l*2 


(30) 


For  the  special  case  where  TiwTj.so  that  zi^Zt^z, 
Eq.(30)  simplifies  to 


V<''i*^o''i''2^^^^v‘* 


(31) 


which  is  precisely  the  same  result  derived  by 
Robinson  in  Ref.3. 


8.  Recursive  Reflection  Transfer  Function 
Relationship 

For  a K-layar  media  systam, states  Uj(t)  and  dj(t) 
have  been  defined  at  the  top  and  at  the  bottom, re- 
spectively, of  the  jth  layer  (sea  Fig.3).l-et  R.(s) 
denote  the  transfer  function  between  Uj  and  ^ at 
tha  jth  Interface;  I . a . , 


i4{u,(t+T.)}  ST.  U,(s) 

■-znrxtrf”"  • 5^  • 


(32) 


When  jwO.wa  obtain  the  reflection  transfer  func- 
tion between  output  y(t)  and  source  m(t);  i.a.. 


i 


sine*  Tg«4. 

W«  shall  now  davtiop  a sinpla  rtcursiva  ralation- 
ship  batwoan  Rj  (s)  and  Rj^^Cs).  Consider  our  ear- 
liar  state  equations  for  u,  and  dj«,i: 

Uj(t-*^j)  - i-jdjCtj+d-rjjuj^jCt)  (3<*) 

'‘j*l^'*^j-rl>  ■ ('♦'j)'lj(«)-r.uj^,(t)  . (35) 

Fron  the  Laplace  transform  of  Eq.(3^).wa  find 

Rj(s)  - r.*(l-rj)Uj.^,(s)/Oj(t)  (36) 

Laplace  transform  Eq.(35)  and  solve  for  Oj(s),to 
show  that 

Oj(s)  - [r.U.^,(s)*,"j^'  0.^,(s)]/(ls.rj)  (37) 

Substitute  Cq.(37)  into  Eq.(36),  and  re-arranqe 
soma  terms  in  the  result inq  expression  to  see  that 

R,(s)  ■ i J (j-K-1,R-2 1,0)  (38) 

which  is  the  desired  result. 

Equation  (38)  can  ba  used  to  compute  the  output  of 
a K-layer  media  system  In  a recursive  manner,  be- 
ginninq  with  a one  layer  system  (I. a., one  layer  on 
cop  of  a basement  layer)  for  which  we  set  jwK-1 . 
we  then  iterate  Eq.(38)  bael(wards,settlnq  JwK-2, 
K-3,...,t,0.  In  order  to  compute  \_^s)  we  need 
R«(s);  but,R,(s)  can  be  obtained  directly  from  the 
very  last  state  equation  In  (9)  ,u,(t*T,)-r,d,(t) , 

IS 

■ ’’k  • (35) 

In  the  special, but  widely  studied  case  of  equal 
travel  t imes .Eqs. (38)  and  (39)  simplify  to 


r ♦t^R  (s) 

R (s)-- L—Jiil— 

J Hrt^R  ,(s) 


j-K-l,K-2 1,0 


• '•k  I 

where  zwe"*"^.  Equation  (i»0)  (or  its  discrete-time 
counterpart, in  which  Laplace  transfer  functions  are 
replaced  by  (-transform  transfer  functions)  is  a 
well-icncMi  result  which  can  be  derived  by  widely 
different  methods  (Refs. 10  and  13,for  example).  Ad- 
ditional ly, these  recursive  relationships  occur  (Ref. 
16)  in  electric  kemal  functions,  maqnetotel luric 
input  impedence  functions, and  electromagnet  I c modi- 
fied kernel  functions. 

That  Eq.(60)  generalizes  to  Eq.(38)  for  non-equal 
t rave  I times  I s be 1 1 eved  to  be  a new  resu 1 1 . 

To  illustrate  the  use  of  Eqs. (38)  and  (39), we  re- 
compute r(s)/N(s)  • R«(s)  for  K>2.  In  that  casa, 

Ri(s)  • ri.and 


^**^2 


'*'’lV2 


1*rg2iR|(s) 


'*'-l'’2Vo'’l*r"oVl*2 

It  Is  much  easier  to  compute  the  reflection  trans- 
fer function  by  the  recursive  reflection  transfer 
function  relationship  of  this  paragraph  than  by  the 
layer  transfer  functions  of  the  preceding  para- 
graph; however,  detailed  information  about  the  up- 
going  and  dawngoing  states  cannot  ba  recovered 
from  the  relationships  in  this  paragraph, whereas 
they  can  be  recovered  from  Eq.(21)  in  the  preced- 
ing paragraph. 

IV.  A RECURSIVE  HALF-STATE  hODEL 

Observe, from  Eq. (38), that  operator  z.^^  appears 
only  as  zf«^  in  Rj(s).  This  suggests  that  a state 
space  modal, which  requires  only  K states, can  be 
developed  from  that  aquation.  We  refer  to  this  mo- 
dal as  a recursive  half-state  model. since  It  is 
obtained  from  the  recursive  reflection  transfer 
function  relationship. 

Theorem  1 . For  the  K-layer  media  system  depicted 
in  Fig. 1 , let 

x_(t)  - col(x,(t),Xj(t) X|j(t))  , (62) 

2 ■ dlag(z,,Z2 2^)  , (63) 

and 

^•(iK-*K-I 

where  e_j  is  the  j-th  unit  vector, and  is  Kx  I.  Then 

TZ‘Vx(t)  - Ax(t)*bm(t)  (65) 

y(t)  - e_*x(t)*rgm(t)  (66) 

where 

b » col(r|^,r|^_,,...,r|)  (67) 

e - col(0,0,...,0,(l-rj))  (68) 


■'■k''k-I 

*''k'’k-2 

*''k'’k-3  ■■ 

■ ■''k'’| 

■''k''o 

(1-rJ-i) 

■'■k-i'’k-2 

*'’k-i'’k-3‘' 

■■'■k-i''i 

*'’k-i'’o 

0 

(l-rj.2) 

*^2^3" 

"'■k-2''i 

■'k-2''0 

0 

0 

(1-^.3)  •• 

■*''k-3'’i 

■'■k-3''o 

\ 0 

0 

0 

•(l-r|) 

-^"o/ 

The  proof  of  this  theorem,  which  contains  a recur- 
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siv«  varsion  of  Eqs.(4S)  and  (46), Is  given  in  Ap* 
pandix  S.  ■ 

The  lafC'iiand  side  of  Eq.(4S)  is  a rsva''se  doubla- 
dalay  operator;  I. a., 

TZ  ^Tx(t)-coI  (x,  (£*27^)  ,X2(t*2T|^.j ) X|^(t*2T,)). 

(50) 

Addi c ional ly, matrix  T Is  a special  permutation  ma- 
trix,with  the  property  T“'"T;  hence, 'q. (4$)  can 
also  be  wri tten.as 

x(t)  - T2^TAx(t)+TZ^Tbm(t)  . (51) 

In  the  special  case  of  equal  travel  timas,Z~zl ,and 
Eq-(5I)  reduces  to 

x(t)  - Ax(t-2T)-r^(t-2T)  . (52) 

Final ly, whereas  the  recursive  half-state  modal  re- 
quires only  K states, we  at  present  can  give  no 
physical  meaning  to  those  states.  In  Section  V.B, 
however, we  relate  the  ^-states  to  the  upgoing 
states  Uj(t),Uj(t),...7u,^(t). 

V.  EQUIVALENT  hOOELS 

Inthis  section  we  present  altemat i ve,but  equival- 
ent,ful  l-state  and  half-state  models.  The  alterna- 
tive ful I -state  models  are  useful  because  they  have 
much  more  revealing  structures  than  the  L-0  full- 
state  model,  and  one  of  these  models  — the  skew 
modal  — leads  quite  naturally  to  a skew  half-state 
model.  That  half-state  modal, whose  states  are  some 
of  the  upgoing  and  some  of  the  dawngoing  states 
from  the  L-0  ful I -state  model, can  be  related  to  the 
Section  IV  recursive  half-state  mode).  In  this  way, 
we  give  physical  meanings  to  the  states  in  the  re- 
cursive half-state  model. 


A,  Ful l-State  Models 
I.  0-U  Model 

The  0-U  full-state  model  is  obtained  from  the  L-0 
full-state  model  by  reordering  the  latter's  equa- 
tions in  such  a manner  that  all  downgoing  states 
are  grouped  together  and  all  upgoing  states  are 
grouped  together.  Let 

d(t)  - col(d,(t),d2(t) d^(t))  (53) 


and 


^(t)  - eol(u,(t),U2(t) 


(54) 


Equations  (9)  and  (3)  can  be  written, in  partitioned 
form  as: 


2 '^(t)  " A^£(t)*Aj£(t)+3m(t) 


Z '^(t)  « Aj£(t)+A^ji(t) 
y(t)  - h'u,(t)T-rgm(t) 


and 

where 


(l*r,) 

0 

0 


0 

0 

(l*rj) 


0 

0 

0 

(I+Tj) 


(55a) 

(55b) 

(56) 

0 \ 

0 
0 
0 


(l^rn-,)  0 


/ 


(57a) 


> 

■ 

"k-I^ 

(57b) 

■ 

< 

diag(rj  ,7^. . . . ,r|^) 

(57c) 

0 

(l-r 

) 0 0 

0 

0 

0 

O 

0 

0 

0 

0 (l-rj)  ••• 

0 

0 

0 0 

(1-r^-, 

'1 

\o 

0 

0 0 

0 

1 

(57d) 

- col  (lTrg,0,0 0) 

(58) 

« col(l-rg,0,0,...,0) 

(59) 

Matrices  Ai,Ai,A],and  A»  are  KMK,and  matrix  oper- 
ator Z is  defined  in  Eq.(43). 

Equations  (55a) , (55b) , and  (56)  comprise  the  D-U 
model.  Matrices  Ai,A2,A],Ak  are  easily  remembered 
for  this  model.  The  0-U  model  has  found  applica- 
tion in  Mendel's  Bremmer  series  decomposition  of 
output  y(t)  (Ref. 15).  That  decomposition  is  dis- 
cussed in  Section  Vtl. 

2,  Skew  Model 


The  skew  full -state  model  is  also  obtained  from  the 
L-0  full-state  model  by  a reordering  or  equations. 
Let  x,,i(t)  and  i,^t)  both  be  K*  I vectors,  where 

x^j(t)  - col (u^(t) ,d2(t) ,Uj(t) ,...)  (60) 
x^2(t)  - col(dj(t),U2(t),dj(t),...).  (61) 


When  K is  even  the  last  elements  of  x^,^t)  and 

are  d,(t)  and  u,  (t) , respectively ; whereas, 
when  K is  odd,  those  elements  are  u,(t)  and  d, (t). 
Equations  (9)  and  (3)  can  be  expressed  In  terms  of 
x,j(t)  and  x,i(t),as 


2 'xj,  (t)  “ 

(62a) 

and 

2 “ Hx^,  (t)+2r'(t) 

(62b) 

vdiere 

y(t)  • il'i,,  (t)+rgm(t) 

(63) 

K • 

1 2 3 4 ,••• 

1 

0 

0 

...\ 

1*r| 

•■■l 

0 

0 

... 

0 

0 

■•3 

'*"3 

... 

0 

0 

’^3 

•^3 

... 

• 

• 

'■7 

6 


r 


1 


I. 


t 


K - t 2 3 


and  and  ar*  defined  In  Eqs.(58)  and  (59)tre- 
sepet  I va  I yT  Matrices  S and  Hare  both  KxK.and  are 
given  above  at  functions  of  K. 


Equations  (62a) , (62b) .and  (63)  conprlse  the  skew 
Model.  It  It  called  a skew  nodal  because  onlycross 
coupling  tenet  appear  on  the  right-hand  side  of  the 
state  equations;  l.e.,tha  rather  unusual  ordering 
^ of  our  2K  states  In  x,/t)  and  x,j(t)  leads  to  a par- 
titioned state  aquatTon  vlth  bTocks  of  zeros  along 
the  main  diagonal . 

8.  Half-State  Models 
1 .Skew  Modal 


Observe, froa  Eqt.(63)  and  (62)  that  y(t)  depends 
only  on  >^,i(t).and  that  Xja(t)  can  easily  be  elimin- 
ated from  the  state  equations  to  give  the  follow- 
ing half-state  state  aquation  for  )^^(t) : 

x^,(t)  - ZGZI^,(t)+ZGZgm(t)  . (66) 

Equations  (66)  and  (63)  constitute  our  skew  half- 
state  modal.  Half-state  vector  £aa(t)  can  be  com- 
puted from  £q.(62b)  once  Xj^(t)  iTas  been  coteputed. 

It  is  interesting  to  compare  Eq.(66)  and  (51).  Ob- 
serve that  the  right-hand  sides  of  these  state 
equations  involve  matrix  operators  ZGZ  and  TZ^T, 
respectively.  If  a physical  element  existed  for 
implementing  a delay  operator  matrix, then  the  skew 
half-state  model  would  need  two  such  elements, 
whereas  the  recursive  half-state  model  would  only 
need  one  such  element  [for  Z*] . This  suggests  that 
the  recursive  half-state  model  Is  minimal  in  terms 
of  "hardware"  requirements. 

In  the  special  case  of  equal  travel  t Imes .where  Z» 
zl,Eq.(66)  reduces  to 

(t)  » GHx^j  (t-2T)+03m(t-2T)  (67) 


which  is  comparable  to  Eq.(52).  We  see, therefore, 
that  only  In  the  case  of  equal  travel  times  do  both 
the  skew  and  recursive  half-state  models  have  the 
saaw  "hardware"  requireaants, In  that  they  both  re- 
quire a double  delay  element, 2t. 

2.  Upgoing  States 

Observe  from  Eq. (56)  that  y(t)  depends  only  on 
jj(t).  While  it  Is  possible  to  el  Imlnate  ^(t)  from 
Tqs.(S5a)  and  (55b)  to  obtain  a half-state  model 
for  ^(t), that  model  is  not  terribly  useful , due  to 
its  Tpparent  complexity.* 

It  is  useful  to  obtain  a direct  relationship  be- 
tween u,(t)  and  ^^(t).  In  Appendix  C,we  show  that 


*u(t)-Z[A^*Aj(l-ZA,)''zA2]u(t)*ZAj(l-ZA,)''z^(t). 


u(t)  - 

where  L,  and  are  permutation  matrices  which  de- 
pend on  whether  K is  even  or  odd.  These  matrices 
are  also  defined  In  Appendix  C. 

From  Eqs.(68)  and  (66) .one  can  develop  a different 
half-state  equation  for  i£(t);  but, that  equation 
also  does  not  appear  to  be  terribly  useful, so  we 
shell  omit  it  here.  Clearly, we  can  compute  u(t)  by 
first  computing  x^x^^)  Then  using  Eq.(6?). 

It  Is  also  possible  to  relate  u(t)  to  x(t).  Wash- 
burn (Ref. 16)  has  proven  the  fol lowing” 

Theorem  2.  Let 

j - 1 

Y,  - H (1+r,)  (69) 

J t-0  ' 

and 

TO)  - It,,  t a 0 . (70) 

I -0  ' ° 

Then, 

*'  • 

Proof:  The  proof  of  this  theroem  is  inductive  in 
nature  and  algebraically  lengthy.  The  general  idea 
Is  to  show  that  £(t)-coj_(c,  (t)  (t) , . ,Cr(t) ) . 

where  C,_j(t)  ■ u,_, (t+T(K-j) )/y,_, ,sat Isf les  the 
recursive  half-state  equation  (45);  i.e.,to  show 
that  for  £(t)  as  just  defined, T2^*T^(t)“A£(t)* 
^m(t) . By  uniqueness, then, ^(t)>x(t) . Details  can 
be  found  in  Washburn  (Ref. 16).  ~m 

(t  is  useful  to  express  Eq.(7t)  in  operator  nota- 
tion. Let 


and 


r - dlag(Y,^,Y|5., Y,) 


0 0 


0 0 


K-  2 

n z, 
i - 1 ' 


M' 


(72) 


(73) 


\ 0 r,  0 - 

\ 1 0 0 .•• 
then, Eq. (71)  can  be  writt 

x(t)  - r’’ 

Theorem  2 is  a very  useful  result;  for, not  only 
does  it  give  us  another  way  to  compute  ^(t)  [i.e., 
u(t)mN*^rx(t) 1 .but  it  also  provides  us  with  mean- 
Tng  for  the  mathematically  defined  half-state  vec- 
tor ,x(t),  in  terms  of  the  physical  half-state,£(t) . 

There  is  an  interesting  pictorial  description  of 
the  relationships  between  X|*x-)  (t)  and  u,(t)  (j" 
1,2....,K).  It  is  obtained  by  writing  Eq.(71)  as 

Uj(t)  - YjX,^^,.j(t-f(j))  . j-1,2 K (75) 

and  is  shown  in  Fig. 6, for  K^k.  On  that  figure.we 


0 0 

0 of 

:en  at 

Nu(t)  . (74) 
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J - 1 

ai«  th«  fact  chat  Y-*  ^ t,,Mn«ra  t:  Is  tha  ith 
■I  I ■ 0 

transmission  coafficlent  (saa  Tabla  1).  States  xi(t) 
through  x»(t)  are  located  at  tha  surface  of  the  <*- 
layer  madia  system, as  mcII  they  should  be  since 
they  are  associated  with  an  input'output  model . The 
dashed  lines  denote  a fictitious  ray  path  linking 
U](c);  and, of  course, we  know  where 
tha  Uj(t)  states  are  located, since  they  are  physl* 
cal  states. The  relationship  between  Uj  and  x,«i_j 
is  one  from  ray  theory.  If, for  example,  X2(t)  Is 
applied  at  the  surface  than, foi lowing  the  direct 
transmission  path  between  X2(t)  and  ut(t),lt  fol* 
lows, from  ray  theory, that  ui(t)“t jtit2X2(t-Ti-T,) . 

In  retrospect,F!g.b  provides  a heuristic  proof  of 
Theorem  2. 


3.  Oowngoinq  States 

Suppose  we  desire  information  about  downgoing 
states  ^(t).  It  is  straightforward  to  show, that 

d(t)  - • (76) 

A procedure  for  computing  £(t)  Is:  (I)  compute 
x,i(t)  using  Eq.(66);  (2)  compute  x^atc)  ^''°m  Eq. 
Te2b) ; and  (J)  compute  £(t)  from  Eq.(76). 


VI . COWmATIQN 

Our  full -state  and  ha  If -state  models  are  causal 
functional  equations,  a class  of  equations  about 
which  we  have  not  been  able  to  find  very  much 
either  in  the  system's  or  mathematics's  literature. 
In  this  section  we  give  some  brief  preliminary 
thoughts  on  simulation  of  causal  functional  equa- 
tions. 

A.  Ful I -State  Model 

Our  discussions  in  this  paragraph  are  directed  at 
simulation  of  the  L-0  full-state  model.  They  are 
applicable,  with  modifications,  to  the  D-U  and  skew 
ful I-state  models. 

A brute  force  approach  to  digital  simulation  is  to: 
discretize  time  (our  independent  variable),  map 
each  Tj  in  some  unique  manner  into  a quantum  zone 
along  the  discretized  time-axis,  and  than  solve  the 
resulting  discrete-time  system  by  means  of  finite- 
difference  equation  techniques.  Appreciable  errors 
may  be  Introduced  in  the  Tj-mapping  step,  unless  a 
very  fine  quantization  Is  used.  To-date,  we  have 
not  tried  this  approach;  but, we  have  not  written  It 
off  either, since  it  is  quite  simple  to  implement. 

A second  approach, which  we  have  begun  to  study  In 
some  detail, is  one  in  which  the  system's  Impulse 
response  is  computed  by*  a branching  proccss/tabla 
look-up  procedure.  Output, y(t) , is  than  computed  via 
convolution  between  m(t)  and  h(t)  [sea  Section  II]. 
Tha  branching  procass/tabla  look-up  procedure  is 
based  on  tha  observation  that  tha  basic  operations 
required  to  compute  the  Impulse  response  are  shift- 
ing and  adding  of  two  non-uniform  saquances  of  im- 
*This  procedure  has  been  developed  by  Nr.  NIchael 
Chan, a graduate  student  In  the  Electrical  En- 
gineering Department  at  the  University  of  Southern 
Cal  I fomia,and  is  reported  on  In  more  detail  In 
Ref.  17. 


pulses.  An  alternate  way  of  viewing  Eqs.(3)  and  (9) 
is  to  consider  what  happens  to  d.(t)  and  uj(t).  A 
careful  study  of  these  equations  leads  to  the  fol- 
lowing transformation  rules: 


<ij.(t) 


rgdo(t)  * y(t) 

r^djCt)  - Uj(t+Tj) 


j-0 

j-1 ,2,...  ,K 

j»0,1 . . . . ,K-1 

(77a) 


(1-ro)U|(t)  -►  y(t)  j-1 

Uj(t)  (I-rj_,)Uj(t)  - Uj.,(frTj_j)  j-2,3 K 

j-I,2,...,K 

(77b) 


We  search  along  the  time  axis  for  a time  at  which 
an  event  (i.e.,an  impulse)  has  occurred.  At  that 
time  point, we  map  all  dj  and  Uj  states  according  to 
Eds.(77a)  and  (77b).  Since  the  right-hand  sides  of 
these  equations  involve  two  time  shifts,  a single 
event  branches  into  two  events.  We  proceed 
along  the  time  axis  looking  up  values  of  dj  and  u. 
at  event  points,  until  we  have  covered  the  domain' 
of  interest.  To  eliminate  costly  computation, we  set 
a lower  bound  on  state  ampl i tudes ,below  which  we 
assume  it  to  be  zero.  Add i t ional ly, i f two  events 
occur  within  a prespecified  tolerance,wve  combine 
results  for  those  events.  The  errors  introduced  by 
these  approximations, as  well  as  storage  require- 
ments are  currently  under  investigation. 


8.  Half-State  Nodels 


Based  on  very  preliminary  results, It  appears  that 
the  recursive  half-state  model  is  computationally 
more  attractive  than  the  skew  half-state  model. 
Coding  the  skew  model  is  difficult  because  we  must 
expand  ZGZH  and  ZGZg^  in  Eq.(66).  We  do  not  know 
general  formulas  for  these  complicated  operator  ma- 
trices. The  recursive  model.on  theother  hand, can  be 
coded  directly  from  Eqs.(45)  and  (46) .making  use  of 
the  very  simple  nature  of  TZ"*T£(t) ,as  given  by  Eq. 
(so).  Either  of  the  two  approaches .described  above 
in  Paragraph  A can  be  used  to  simulate  the  recur- 
sive model.  We  are  presently  studying  such  ap- 
proaches . 


VII.  CONCLUSIONS 

We  have  developed  state  space  models  for  lossless 
layered  media  which  are  described  by  the  wave  equa- 
tion and  boundary  conditions.  Our  models  are  for 
non-equal  one-way  travel  times, and  are  therefore 
more  general  than  traditional  transfer  function 
models.which  are  usually  for  layers  of  equal  one- 
way travel  times.  Our  state  rpace  models  represent 
a new  class  of  equat Ions .which  we  call  causal  func- 
tional equations.  These  equations  are  I inear, time- 
invariant, continuous-time  equations  with  multiple 
time  delays.  The  Impulse  response  of  our  system  is 
a saquence  of  unequally  spaced  Impulse  functions. 

We  have  developed  full-state  models.which  require 
ZK  states  for  a K-layer  media  system, and, have  also 
developed  half-state  models.which  require  only  K 
states.  Add i t ional I y, we  have  generalized  certain 
transfer  funct ions .which  aopear  in  the  geophysics 
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1 1 Caratura, from  layarad  madia  with  tquai  travel  last  author  wishes  to  thank  the  Hughes  Feilowship 

times  to  layered  media  with  unequal  travel  times.  Foundation  for  his  support. 


Now  that  we  have  developed  state  space  models  for  a 
K-layer  madia  systaai.much  work  remains  before  us. 
Since  we  have  been  led  to  a new  class  of  equations, 
causal  functional  aquations, they  must  be  studied 
not  only  from  a simulation  point  of  view, but  also 
from  a system  theoretic  point  of  view.  Efficient 
coatputat ional  methods  must  be  dove i oped , and  notions 
such  as  obsarvabi 1 1 ty,cantrol labi I ity ,and  identifi* 
ability  must  be  expanded  to  this  new  class  of  equa- 
tions. Work  is  presently  underway  in  these  areas. 

Additionai  areas  of  study, which  are  also  underway, 
all  deal  with  what  can  bo  done  with  the  state  space 
models.  These  studies  include:  (1)  extending  the 
modal  to  include  absorption  and  non-normal  inci- 
dence effects;  (2)  identifying  reflection  coeffici- 
ents and  travel  times  using  a recursive  layer- 
stripping procedure;  and  (3)  developing  minimum- 
variance  state  estimators. 

In  closing, we  wish  to  sumearlza  a very  Interesting 
decomposition  of  autput,y(t) ,which  was  developed  by 
Handel  (kef. IS), but  is  originally  due  to  Braniier 
(Ref. 18).  This  decomposition  was  vary  easy  to  de- 
velop using  our  state  space  models.  The  decomposi- 
tion,which  we  refer  to  as  a canonical  Bremmer 
Series  degompos i t i on , i s : the  complete  output, y(t), 
from  a K-layer  media  system, which  is  comprised  of 
the  superposition  of  primaries,secondaries^terti- 
aries,atc. ,can  be  obtained  from  a single  state 
space  model  of  order  2K  — the  complete  modal  ~ or 
from  an  infinite  number  of  modeis,each  of  order  2K, 
the  output  of  the  first  of  which  Is  just  the  pri- 
maries,the  output  of  the  second  of  which  is  just 
the  secondaries, etc. 

By  thinking  of  y(t)  as  the  superposition  of  Its 
constituents  — primaries, secondaries, etc. , — we 
can  write  it  as 


y(t)  - J y.(t)  . (78) 

j-1  J 

Equation  (78)  is  the  Bremmer  Series  decomposition 
of  y(t).  Bremmer  shows  how  to  compute  the  constitu- 
ents,y ,(t)  .from  integral  equations  which  relate 
yj(t)  to  y^^(t).  Hendel ,on  the  ocher  hand, shows  how 
CO  compute  the  yj(t)  as  depicted  In  Fig. 7.  Input 
m(c)  drives  a state  space  primaries  model .whose  up- 
going  states  drive  a state  space  secondaries  model, 
etc.  The  0-U  full-state  model  Is  most  appropriate 
for  characterizing  Che  Bremmer  Series  decomposition. 
Further  details  on  the  structure  of  the  n-arles  mo- 
del as  well  as  a proof  of  validity  of  the  decompo- 
sition are  given  In  Ref. IS. 
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appendix  a.  Oerivation  of  Laver  Transfer  Functions 

We  direct  our  attention  at  the  upgoing  and  downgo- 
ng  states  in  each  layer  [Eq.(9)].  it  is  straight- 
forward to  show  that  the  complete  set  of  2K  Laplace 
transformed  state  equations  can  be  written, as 


F,X,  - G,M+H,X, 
1—1  1 I-^ 

F.X^  - G,X,rH,X^ 


F X 


'•K-liK-2'*^K-liK 


^lA(  “ ‘'lA.  I 


(A-1) 


where  explicit  dependence  of  all  quantities  on  s 
has  been  omitted  for  notatlonal  simplicity.  Ma- 
trices F,  ,H,,and  G,  are  defined  in  Eqs . (24) , (25) , 
and  (26)  and  (27) .respectively.  From  the  last  equa- 
tion In  (A-I),wa  see  that 


’'k  “iA-i 


^ w*'g  X 
K V)k-1 


(A-2) 


Substitute  this  equation  into  the  next  to  the  last 
aquation  In  (A-l),to  show  that 


^-1 


K'  ^K-I^-2 


(A-3) 


which  can  also  be  written, as 


where 


■ -kVk-A-2 


U ^ F — H 
K-I  K-l  K-l  K K ■ 


(A-4) 

(A-S) 


Proceeding  in  a similar  manner, it  is  easily  shown 
that  Eq.(A-l)  can  be  expressed  In  terms  of  matrix 
W)  ,as 


X,  - W,  G,M 

A ■ “2 '“2^.1 


(A-6) 


where  Is  defined  in  Eq.(23). 

In  order  to  obtain  the  state  transfer  function  in 
Eq.  (21 ) ,subst i tute  X_i  lnto£i,£j  into  Xj,...,  end 
X^.ilnto  Xj,.  In  orde’r  to  obtain  the  reflect  ion 
transfer  function  Eg. (22) .first  express  Y(s),from 
Eq.  (3) , in  te-ms  of  )^t(s),as 

Y(s)  • (l-hg)  (1 ,0)X^I  (s)+rgM(s)  , (A-7) 


and  then  substitute  Xi,  from  Eq. (A-6) , Into  Eq.(A-7). 


APPENDIX  B.  Oerivation  of  Recursive  Hal f-State  Model 

For  the  purposes  of  deriving  the  recursive  half- 
state  model, it  Is  convenient  to  use  a recursive 
transfer  function  comparable  to  Eq.(38)  which  iter- 
ates in  a forward  direction  instead  of  a backwards 
direction,  in  this  way, a one-layer  system  will  be 
associated  with  subscripted  "one"  quantities  rather 
than  subscripted  "K"  quant i t ies .etc.  Consider  Eq. 
(38)  for  j«0: 
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R^(») 


(B-I) 


Bcnuoibar  cha  layars  of  our  K-layar  madia  syston  as 
dapictad  In  fig. 3b, in  which  tha  top  layar  is  now 
tha  Kth  layar._haking  tha  transformations  rj-r,  , 

* *«  .3«  (*)  ■*■  3«(»)  .*nd  Ri  (s)  - 3,_i(s)  .Eg.  (B-1)  can 
ba  wri ttan  as 


ys) 


(8-2) 


Rafleccion  transfar  function  !'<($)  can  ba  obtainad 
in  a racursiva  mannar  from 


R.(s) 

J 


^ S-i  ' . J-1.2.. 

t.r.z.Rj.,(s) 


...K 


(B-3) 


Wa  prova  Theorem  1 by  first  stating  and  proving  the 
truth  of  a recursive  half-state  modal  for  tha 
laver-ravarsed  system  in  Pig. 3b, and  then  transform- 
ing  variables  from  that  modal  to  tha  Pig. 3a  situa- 
t ion. 


Theorem  3 • Por  tha  K-layar  media  system  depleted 
In  Pig. 3b, I ec 


x-'(t)  - coKx  (t),x,(t),...,x  (t))  (8-l») 

and  ~ ■* 

Z.  - diag(lj,22 Z.)  (B-5) 

where  j«l,2,...,K.  Then, 

Zj^x;*(t)  ■ A-^x;*  (t)+b_'*mj  (t)  (B-6) 

y*(t)  • (t)+r^.m^.  (t)  (B-7) 


where  m.(t)  and  y'i(t)  denote  the  input  and  output, 
respectively, at  the  surface  of  a j-iayer  media 
system, and 


1 

c 


-'•o''i 


(1-rJ) 


(B-8) 


c_-^  - col(0.0,....0,(l-r^))  . 

j-' 


(B-9) 


(B-10) 


(B-n) 


Proof ; Our  proof, which  is  by  induct  ion, is  based 
on  the  block  diagram  decomposition  of  Eg. (B-3), de- 
picted in  Fig. 9.  To  begin, we  develop  the  initiali- 
zation in  Eg. (B-8);  than, we  demonstrate  the  truth 
of  tha  theorem  for  j’‘2;  and, f .nal  ly, assuming  the 
truth  of  tha  theorem  for  j-l,wa  demonstrate  its 
truth  for  j . Because  the  steps  for  j“2  and  j are 
so  similar  in  nature, we  omit  the  j^Z  proof. 

(1)  Initializations  (jwi) 

From  F1g.9,wo  find, for  j»1,that  [R^o(s)  • r,] 

z,^x'(t)  • -rgrjx\c)+rgmj  (t)  (8-12) 

and 

y*(t)  ■ (1-rj)x’(t)*rjm^(t)  . (8-13)- 

Comparing  Egs . (8-1Z)' and  (B-i3)  with  Egs.(8-6)  and 
(8-7) .respectively, it  is  clear  that  A',b‘,and  c' 
are  as  defined  In  Eg. (B-8). 

(2)  Inductive  Step 

Observe  that  the  model  of  the  (j-l)-Iayer  media 
system  Is  embedded  In  the  forward  path  of  the 
feedback  loop  in  the  model  for  the  j-layer  media 
system.  We  hav^  labeled  the  input  and  output 
points  of  the  R]_i(s)  block  as  h,,^(s)  and  Vil-i(s), 
respectively.  Observe, that 

2:V(t)  - yJ''(t)  ; (B-lk) 

but, by  assumption, yi-Kt)  satisfies  Eg.  (B-7);  hence 
Zj^x-*(t)  • c'"**'*'^*' (t)+r^._,m^._j(t)  . (B-15) 

Prom  Pig .9, however 

mj_^(t)  ■ m^ (t)-rjX^ (t)  (8-16) 

whereupon  Eg. (B-15)  becomes 

2j^x'*(t)  - £‘-‘  ’x-*  ’(t)-r^.rj_|X-'(t)wrj_jm^.(t) 

(B-17) 

which  is  our  state  eguation  for  x^ . 

Next,  we  must  combine  Eg. (8-17)  with  our  assumed 
state  eguation  for  x^“\Eg. (B-6) ; but , we  must  re- 
place m<.i(t)  in  the  latter  eguation  by  £g.(B-16). 
The  eguation  for  x;"^  is 

2j?,£'*’(t)  - Aj*'xj*’(t)-rj.b-''’x-'(t)+b-''’m^.(t)  . 

■“  ■'  (B-18) 

Clearly  now,Egs . (B-l8)  and  (B-17)  can  be  combined 
and  expressed  as  in  Eg. (B-6) , where  A^  and  ^ are 
given  by  Egs. (8-9)  and  (B-IOKrespect ively. 

Finally, from  Pig. 6, we  see  that 

yj(t)  - (1-rj)x-'(t)+r.m.(t)  (B-I9) 

and, this  can  be  expressed  as  in  Eg. (B-7) .where  £* 
is  defined  in  Eg.(B-ll).  ■ 

While  Egs.(B-9)and  (8-10)  are  interesting  in  their 
own  right, they  do  not  reveal  the  intrinsic  detai  led 
structures  of  A^  and  b^* . It  is  a straightforward 
matter  to  iterate  these  two  squat  ions, using  their 
starting  values  from  Eg. (B-8), to  show  that 
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• '“'('•o-'". 


It  i>  straightforward  to  show  that 

!!,<•'  • "A,'" 

uJjCt)  - L^2(t) 


■'■o''2 

-^0^3 

-Vj-1  -Vj  \ 

where 

L,  and 

L»  are 

both  K/2  » 

: K matrices 

,and 

-•/z 

. 7. 

•^'•3 

and 

L 

e 

• (s.,-o 

.2^.0.... 

•*K/2'-’ 

(C-8) 

{l-r‘) 

■^2^3 

S, 

* (£*•) 

,0,ej,... 

’-’2</2^  • 

(C-9) 

0 

(1-rj) 

’Vj-l  ■"3''j 

For  (C 

odd , 

^ y 

ij/t) 

^col  (u^l 

[t),Uj(t), 

....u^t) 

(dj(t),d^(t),, 

....d,.,(t)) 

0 

0 

•••  ^'---j-lj  *"j-1"j/ 

-col (u* 

|(t)(d,*, 

(t)) 

(C-IO) 

Proof  of  Thaorem  1 ! From  Ftg.S.wa  saa  that  tha 
state  space  modal  for  tha  desired  system, Fig. 8a, 
can  be  obtained  from  the  state  space  modal  for  tha 
lavar-revarsed  system,FI g .3b,by  means  of  tha  fol* 
lowing  transformat  Ion  of  variables: 

*l-*Vl  + l ''  • 

and  ■*  ■' 

•■i  * '"v.i  • (B-23) 


In  order  to  obtain  A,b,and  c In  Eqs. (49) , (47) ,and 
(48) , respect lvely,appTy  Eq.r8-23)  to  Eqi.(B*2l), 
(a-20),and  (8- 1 1 ) ,and , then  sat  J“K.  In  order  to 
obtain  the  laft-hand  side  of  Eq.  (4;)  ,apply  Eq.(8-22) 
to  2*,*  (In  Eq.  (8-6)1  •”<1  It  Is, of 

course  understood  that  quantities  In  Theorem  3 
which  are  superscripted  or  subscripted  "K"  [l.a., 
^‘(t)]  are  tha  same  as  those  unsuparscr Ipted  or  un- 
subscrlptad  quantities  In  Theoram  1. 

APPENDIX  C.  Derivation  of  Relationships  between 

a , (t)  and  u(t) 

— s I — 

Hare, we  shall  show  that  there  exist  permutation 
matrices  L,  and  L«, whose  dimensions  and  structures 
depend  on  whether  K Is  even  or  odd, such  that 

u(t)  . L;LaX„(t)*L'L^2(‘)  ■ (=*'> 

Substitute  Eq. (62b), for  ^,e(t),lnto  Eq.(C-l)  to 
show  that 

u(t)  - (L;i.^+l.;L^ZH)x^,(t)+^L^Zam(t)  . (C-2) 

It  Is  straightforward  to  show, by  direct  substitu- 
tion of  tha  defining  equations  for  L,,Z,and  ait6at 
L^LtZa  * 0;  hence, 

u(t)  • (L;L^+I.^L^ZH)x^,(t)  (C-3) 

which  Is  tha  desired  result. 

Now  wo  return  to  tha  decomposition  of  u_(t)  'h  Eq. 
(C-I).  Let  xji(t)  and  denote  rear~angad  x,i(t) 

and  vectors.  For  K even. 

(t)licol(u^  (t),Uj(t). . . .01^.,  (t)  |dj(t), . . . .d|^(t)) 

•coKu*  (t)l^,  (t))  ^ (C-4) 

and 

xj2(t)ieol  (uj(t),Uj^(t), . . . ,U|^(t)  |dj  (t) du_|  (t) ) 

-col  (u*,(t)  |^,(t))  . (C-5) 


-col  (u*2(t)  l^j(t))  . 


xJ^t)^col(oj(c),u^(t), . ..  ,u^.,(t)  |dj(t),dj(t) 

-col  (u*2(t)  |^2(t))  . (C-ll) 

Equations  (C-6)  and  (C-7)  are  applicable  In  this 
case  also;  but, now  L,  Is  (—^  * K and  is  * 

K,and 




S ■ 2«-r-*  • 


In  Eqs.(C“8)  and  (C-9)  .Cj  are  yx  i unit  v«ctor5 ; 

but, In  Eq.(C-l2),  ej  are  I, and, In  Eq.(C-I3). 

K-l  , 

Oj  are  -J— * I • 

Pagardlass  of  whether  K Is  even  or  odd,  one  can 

shcM  that 

u(t)  - 4^,  (t)+L|^uJ,  (t)  . (C-I4) 

Substitute  Eqs. (C-6)  and  (C-7)  Into  Eq.(C-l4)  to 

Obtain  tha  assumed  decomposition  for  u(t)  In  Eq. 

(C-I). 
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Figure  I.  System  of  K Layered  Media 


Figure  3. 


Reflected  and  Transmitted  Waves 


Interface.k.  From  Eq.(8), 

d4(t-T,)  ^d,(t). 


at 


Figure  4.  Summary  of  Past  (top  path)  and  Present 
(lower  path)  Work.  Dashed  Block  Denotes 
Work  too  Preliminary  to  be  Reported  on 
in  this  Paper. 
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Raflected  and  Transmitted  Uavas 
at  (a)  Surface  and  (b)  Intar 
face  K. 


System  of  K layarod  madia: 

(a)  first  layer  Is  Layer  1 and 
tha  last  layer  Is  Layer  K; 

(b)  first  layer  Is  Layer  K and 
tha  last  layer  Is  Layer  1. 


Relationship  between  j((t)  and  u(t). 
Updoing  states  are  located  (circles) 
at  top  of  each  layer.  Dashed  I Ines 
denote  transmission  paths  for  x 
states.  States  are  separated  along 
the  horizontal  axis  for  purposes  of 
clari ty . 


B^lock  Diagram  Representation  of 
Rj(s)  fora  J-layer  system, j>1 ,2 

...,K.  For  jw1 , disregard  the 
M.(s)  and  T*^(s)  labels. 
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Canonical  Bremmer  Series  decompo- 
sition of  a seismogram  signal, y(t) 
Vector  u denotes  the  collection 

of  K upgoing  states  from  the  n- 
arles  model . 
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